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Abstract 

We explore the singlet scalar dark matter (DM) from direct detections and high energy neutrino signals 
generated by the solar DM annihilation. Two singlet scalar DM models are discussed, one is the real singlet 
scalar DM model as the simple extension of the standard model (SSDM-SM) with a discrete Z2 symmetry, 
and another is the complex singlet scalar DM model as the simple extension of the left-right symmetric 
two Higgs bidoublet model (SSDM-2HBDM) with P and CP symmetries. To derive the Sun capture rate, 
we consider the uncertainties in the hadronic matrix elements and calculate the spin-independent DM- 
nucleon elastic scattering cross section. We find that the predicted neutrino induced upgoing muon fluxes 
in the region 3.7 GeV < mi, < 4.2 GeV slightly exceed the Super-Kamiokande limit in the SSDM-SM. 
However, this exceeded region can be excluded by the current DM direct detection experiments. For the 
SSDM-2HBDM, one may adjust the Yukawa couplings to avoid the direct detection limits and enhance 
the predicted muon fluxes. For the allowed parameter space of the SSDM-SM and SSDM-2HBDM, the 
produced muon fluxes in the Super-Kamiokande and muon event rates in the IceCube aie less than the 
experiment upper bound and atmosphere background, respectively. 
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I. INTRODUCTION 

The existence of dark matter (DM) is by now well confirmed . The recent cosmologi- 
cal observations have helped to establish the concordance cosmological model where the present 
Universe consists of about 73% dark energy, 23% dark matter and 4% atoms Js]. Understanding 
the nature of dark matter is one of the most challenging problems in particle physics and cosmol- 
ogy. Currently, many DM search experiments are under way. These experiments can be classified 
as the direct DM searches and the indirect DM searches. The direct DM detection experiments 
may observe the elastic scattering of DM particles with nuclei. The indirect DM searches are 
designed to detect the DM annihilation productions, which include neutrinos, gamma rays, elec- 
trons, positrons, protons and antiprotons. In addition, the collider DM searches at CERN LHC are 
complementary to the direct and indirect DM detection experiments. 

The indirect DM searches are usually independent of the direct DM searches. Namely, one can 
calculate the DM annihilation signals when the thermal- average of the annihilation cross section 
times the relative velocity (crv) and the DM annihilation productions are known. It is worthwhile 
to stress that the DM annihilation signals from the Sun (or Earth) depend on both the direct DM 
detection and the indirect DM detection. When the DM particles elastically scatter with nuclei in 
the Sun, they may lose most of their energy and are trapped by the Sun yj]. The solar DM capture 
rate is related to the DM-nucleon elastic scattering cross section. These trapped DM particles 
will be accumulated in the core of the Sun due to repeated scatters and the gravity potential. 
Therefore the Sun is a very interesting place for us to search the DM annihilation signals [|4j-l9|]. 
The DM annihilation rate in the Sun depends on (crv) and the solar DM distribution. If the DM 
annihilation rate reaches equilibrium with the DM capture rate, the solar DM annihilation rate 
only depend on the DM-nucleon elastic scattering cross section. Due to the interactions of the DM 
annihilation products in the Sun, only the neutrino can escape from the Sun and reach the Earth. 
These high energy neutrinos interact with the Earth rock or ice to produce upgoing muons which 
may be detected by the water Cherenkov detector Super-Kamiokande (SK) jlOl] and the neutrino 
telescope IceCube jlll] . 

In this paper, we explore the singlet scalar dark matter from direct detections and high energy 
neutrino signals via the solar DM annihilation in two singlet scalar DM models. One is the real 



singlet scalar DM model as the simple extension of the standard model (SSDM-SM) 1112 - 



m and 



another is the complex singlet scalar DM model as a simple extension of the left-right symmetric 



2 



two Higgs bidoublet model (SSDM-2HBDM) il7l4l9ll. In the SSDM-SM, a real singlet scalar 



S with a Zj symmetry is introduced to extend the standard model. Although this model is very 



simple, it is phenomenologically interesting [Il2l-ll6ll. In the SSDM-2HBDM, the imaginary part 



S D of a complex singlet scalar field S = (S + iS d)/ V2 with P and CP symmetries can be the 
DM candidate [jiol . The stability of So is ensured by the fundamental symmetries P and CP of 



quantum field theory. In Refs. [|16l1 and lioll . we have calculated the spin-independent DM elastic 
scattering cross section on a nucleon. In fact, one should consider the uncertainties in the DM 
direct detection induced by the uncertainties in the hadronic matrix elements. Here we consider 
these uncertainties and recalculate the spin-independent DM-nucleon elastic scattering cross sec- 
tion. Then we calculate the neutrino fluxes from the singlet scalar DM annihilation in the Sun 
and the neutrino induced upgoing muon fluxes in the Super-Kamiokande and IceCube. This paper 
is organized as follows: In Sec. II, we outline the main features of the SSDM-SM and SSDM- 
2HBDM, and give the DM-nucleon elastic scattering cross section. In Sec. Ill, we numerically 
calculate the diff'erential neutrino energy spectrum generated by per DM pair annihilation, the DM 
annihilation rate in the Sun and the neutrino induced upgoing muon fluxes. Some discussions and 
conclusions are given in Sec. IV. 

n. CONSTRAINT ON SINGLET SCALAR DARK MATTER FROM DIRECT DETECTIONS 

A. The real singlet scalar dark matter model as an extension of the SM 

In the SSDM-SM, the Lagrangian reads 

X = XsM + \d^Sd^S -^S^-^S'- AS^H^H , (1) 

where H is the SM Higgs doublet. The linear and cubic terms of the scalar S are forbidden by 
the Z2 symmetry 5^-5. Then S has a vanishing vacuum expectation value (VEV) (5 ) = 
which ensures the DM candidate S stable. As describes the DM self-interaction strength which 
is independent of the DM annihilation. After the spontaneous symmetry breaking (SSB), one can 
obtain the DM mass = ml + A v^^ with vew = 246 GeV. The SSDM-SM is very simple and 
has only three free parameters: the DM mass m^,, the Higgs mass m/, and the coupling A. In terms 
of the observed DM abundance Qdm^^ = 0.1123 + 0.0035 [3|], one can calculate the coupling A 
for the given mo and m/,. Here we take m/, = 125 GeV n2(M and 1 GeV < m^ < 200 GeV for 



illustration. As shown in Fig. [T] the observed DM abundance requires A ~ (9(10 ^ - 1). It is well 



known that the annihilation cross section a will become larger for the same coupling when the 
annihilation process nears a resonance. This feature implies that there is a very small coupling 
when 0.8 m/, < Imjj < nih as shown in Fig. [U This region is named as the resonance region in the 
following parts of this paper. 
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FIG. 1: The predicted couplings A in the SSDM-SM (left panel) and Ai^d in the SSDM-2HBDM (right 
panel) as a function of the DM mass nio from the observed DM abundance Odm^^ = 0.1123 ± 0.0035 |13|]. 
The vertical dashed line with arrowhead in the left panel shows the excluded region from the potential's 
global minimum, perturbativity and DM relic density. 



Using the predicted A from the observed DM abundance, one can calculate the spin-independent 
DM-nucleon elastic scattering cross section [12 ill 
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where m„ is the nucleon mass and / = (7/9) 'Z,q=u,d,s fj-q + 2/9. In terms of the relevant formulas 
in Ref. [SI], one can calculate the parameters fj^ and obtain / « 0.56 + 0.17. On the other hand, 
the lattice results imply / « 0.29 + 0.03 where we take the strange-quark sigma term 16 MeV < 
cr, < 69 MeV ||22|| . Therefore we adopt 0.26 < f < 0.73 for the following analyses. The 
authors in Re f. [il5i1 have discussed that the light DM particle S can explain the DAMA ll23n 
and CoGeNT ll24ll experiments. Here we consider the latest experiment limits and recalculate the 
spin-independent DM-nucleon elastic scattering cross section crf-,^ with 0.26 < f < 0.73. Notice 
that crf^ is not sensitive to the Higgs mass in the low DM mass range. As shown in Fig. [2] (top 
left panel), the predicted erf in the region 6 GeV < m/) < 8 GeV and / > 0.60 well fit the 



common region of the DAMA and CoGeNT ll25ll . However, the recent CDMS II ll26ll disfavors 



the CoGeNT+DAMA region. We find that the CDMS II Um, CDMS (shallow-site data) MM, 
CRESST Urn and TEXONO m\ can exclude the / > 0.63 region for 1 GeV < mo < 10 GeV. 
The / = 0.63 case has been shown as the blue solid line in the top left panel of Fig. [2l The latest 
XENON 100 im may exclude 7 GeV < mo < 52 GeV and a narrow region 65 GeV < fnn< 80 
GeV even if we take / = 0.26 as shown in Fig. |21 The future experiments CDMS 100 kg plQ and 
XENONIT 1I32I1 can cover most parts of the allowed parameter space. As the DM mass increases, 
new DM annihilation channels will be open which means that the predicted A from the DM relic 
density will quickly decrease. Therefore a kink around the bottom quark mass ^ mi, = 4.2 
GeV occurs in the top left panel of Fig. |2] When the DM mass approaches the resonant point 
= 62.5 GeV for m/, =125 GeV, one can obtain a very large thermally averaged annihilation cross 
section (crv) for the given A. In order to derive the correct DM relic density, we have to require 
that the coupling A is very small. In this case, a very small crfj^ around the resonant point can be 
obtained. 

The SSDM-SM also suffers other constraints except for the direct detections, such as the po- 
tential's global minimum at (h) = vew and (5) = requires \A\ < ^As /IniiJvEw + "^d/^ew alM- 



Since the perturbativity implies 6As < An, one can derive \A\ < ^Jn/3mh/vEw + "^o/^ew Then 
we find the desired DM relic density can exclude < 4.1 GeV for m/, = 125 GeV. The vertical 
dashed line with arrowhead in the left panels of Figs. [T] and |2] shows the excluded region. In Ref. 
ll^ . the authors have given the lower bounds on for several typical As based on the one-loop 
vacuum stability and the observed DM relic density. 



B, The complex singlet scalar dark matter model as an extension of the 2HBDM 



We begin with a brief review of the 2HBDM described in Ref. [|l7|l. The model is based on the 
gauge group SU(2)l ® S U(2)r ® U(1)b-l- The left- and right-handed fermions belong to SU{2)l 
and SU(2)r doublets, respectively. The Higgs sector contains two Higgs bidoublets (p (2,2*,0),^ 
(2,2* ,0) and a left(right)-handed Higgs triplet A^^) (3(1),1(3),2) with the following flavor contents 
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The introduction of Higgs bidoublets <p and;^' can account for the electroweak symmetry breaking 
and overcome the fine-tuning problem in generating the spontaneous CP violation in the left-right 
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FIG. 2: The predicted DM-nucleon elastic scattering cross section af^ in the SSDM-SM (top row) and 
SSDM-2HBDM (bottom row) for 1 GeV < mo < 200 GeV. The dashed Unes indicate the current experi- 
mental upper bounds. The short dotted lines in the right panels denote the future experimental upper bounds 



from the CDMS 100 kg and XENONIT tM- The blue solid line in the top left panel describes the 



/ = 0.63 case. The black hatched region corresponds to a combination of the DAMA and CoGeNT 12511 . 
The vertical dashed line with arrowhead in the top left panel shows the excluded region from the potential's 
global minimum, perturbativity and DM relic density. 



symmetric one Higgs bidoublet model. Meanwhile it also relaxes the severe low energy phe- 
nomenological constraints nlln . Motivated by the spontaneous P and CP violations, we require P 
and CP invariance of the Lagrangian, which strongly restricts the structure of the Higgs potential. 
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The most general potential containing only the (p and Al,r fields is given by 

%A = -jU?Tr(0V) - )"2[Tr(0V) + Tr(^0t)] - //^[Tr(ALA[) + Tr(ARA^)] 
+ii[Tr(<^V)]' + ^2{[Tr(^V)]' + [Tr(^0^)]'} + i3[Tr(^V)Tr(^0^)] 
+i4{Tr(0V)[Tr(0V) + Tr(^0^)]} 

+Pi{[Tr(ALA[)]2 + [Tr(ARA;)]2} +p2[Tr(ALAL)Tr(A[A[) + Tr(ARAR)Tr(A^A^)] 
+P3[Tr(ALA[)Tr(ARA;)] +p4[Tr(ALAL)Tr(A;A^) + Tr(A[A[)Tr(ARAR)] 
+aiTr(0V)[Tr(ALA[) + TrCARA^)] + a2Tr[(0V) + (^0^)]Tr[(ALA[) + (ArA^)] 
+Q'3[Tr(<^0"fALA[) + TrC^^ARA^)] 

+;8i[Tr(0AR0l^A[) + Tr(</.tAL<^A;)] +;e2[Tr(MR0U[) + Tr(0l^ALMR)] 
+;83[Tr(<^AR0^A[) + Tr(0tALMi)], (4) 

where the coefficients /I,, p,, a, and yS, in the potential are all real as all the terms are self- 
Hermitian. The Higgs potential "V^a involving^ field can be obtained by the replacement^ <r^ (p 
in Eq. dH). The mixing term 'Vx<i>s can be obtained by replacing one of by ;^ in all the possible 
ways in Eq. dH). After the SSB, the Higgs multiplets obtain nonzero VEVs 



where /Ci, K2, Wi, W25 and are in general complex, and k = a/I/CiP + I/C2P + |wiP + |w2p « 246 
GeV represents the electroweak symmetry breaking scale. The value of vr sets the scale of left- 
right symmetry breaking which is directly linked to the right-handed gauge boson masses, is 
subjected to strong constraints from the K, B meson mixing as well as low energy electroweak 
interactions. The kaon mass difference and the indirect CP violation quantity set a bound for 



Vr around 10 TeV [|33|]- In general, the 2HBDM includes three light neutral Higgs bosons and a 
pair of charged light Higgs particles, whose masses are order of the electroweak energy scale. For 
simplicity, we consider /<:2 ~ W2 ~ 0. Then one can derive three light neutral Higgs bosons: h,H,A 
from (/ij and;^^j, and a pair of charged light Higgs particles: h- from;^^*. For a concrete numerical 
illustration, we choose all the masses mh, m^, mh± = 180 GeV and m/, = 125 GeV. 

In the 2HBDM, the P and CP symmetries have been required to be exactly conserved before 
the SSB, thus the discrete symmetries P and CP can be used to stabilize the DM candidate. In 
the framework of 2HBDM with a complex singlet scalar S = {S^ + iSi))! ^ (SSDM-2HBDM), 



we have considered this possibility in Ref. lll9ll . The P and CP transformation properties of the 
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TABLE I: The P and CP transformation properties of the Higgs particles and their gauge-invariant combi- 
nations. The and "-" denote even and odd, respectively. 

Higgs particles and their gauge-invariant combinations have been shown in Table HI It is clear that 
the odd powers of (S - S*) are forbidden by the P and CP symmetries. This hidden discrete Z2 
symmetry on 5^ is induced from the original P and CP symmetries. With the help of this hidden 
Z2 symmetry, one may derive (S^) = or (S o) for the VEV of S d- Since the (So) + 
case means that S d may decay and can not be the DM candidate, we require that S obtains a real 
VEV {S ) = v^l V2. Although both P and CP are broken after the SSB, there is still a residual Z2 
symmetry on S d- Therefore 5/) is a stable particle and can be the DM candidate. We have checked 
that the P and CP transformation rules for S defined in Table H] is actually the only possible way 
for the implementation of the DM candidate. 

For the annihilation cross section of approximately weak strength, we expect that the DM mass 
is in the range of a few GeV and a few hundred GeV However, the mass mo of is related to the 
LR symmetry breaking scale ~ 10 TeV To have a possible light DM mass, we may consider 
an approximate global U(l) symmetry on 5, i.e. 5 e'^S . Then the P and CP invariant Higgs 
potential involving the singlet S is given by 

7 2 

= -As S * + An(S S*f + Y, \dS 5 - ^(5 - 5 *)^ (6) 

1=1 

where O, = Tr(0V), O2 = Tr(0t0 + o, = Trix'^x), O4 = Trix'^x O, = Tr(0V +xH\ 

Oe = Tr(x'^^ + <j)^x) and O7 = Tr(A[AL -I- A^A/?). Only the last term explicitly violates U{1) 
symmetry. After the SSB, S obtains a real VEV Va-I V2. Then one can straightly derive 

% = -fViSl + Iv^S^ + Slf- v',] + -f^Sl + Iv^S^ + vl + Sl){Oi - (0,)) + , (7) 

(=1 

where we have used the minimization condition //^ = /Idv^ + JjiAt^oiOi) from the singlet S a- to 
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eliminate the parameter hd- The terms proportional to odd powers of 5 d are absent in Eq. d?]) 
which implies S d can only be produced by pairs. Notice that the mass term oi S d should be 
absent with an exact global U{1) symmetry. As discussed in Ref. [|l9ll . the explicit breaking of 
this U{1) symmetry can explain the naturalness of a light DM mass m^, but it does not destroy the 
stability of the DM candidate S o- For the VEV of So-, we require Vo- > ~ 10 TeV » k which 
means the mixing angles between So- and other neutral Higgs bosons in the SSDM-2HBDM are 
small and the mass of 5 o- is very heavy. 
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FIG. 3: Feynman diagrams for the DM annihilation in the SSDM-2HBDM. 

For the DM mass, we take 1 GeV < mo < 200 GeV. In this case, the possible DM anni- 
hilation products are //, W\W\IZ\Zi, Wfh^/Zi(h,H,A), h^h' and any two of the three neutral 
states (h,H,A) as shown in Fig. O Here Wi and Zi denote the SM gauge bosons. For cubic and 
quartic scalar vertexes, we assume they are the same as those in the one Higgs bidoublet case 



II19I1 . Namely, the vertexes of S dS D(h,H,A) and S dS oih, H,A/h^){h, H,A/h ) are set equal to 



-z/li £)Vew and -iAi^o, respectively. Similarly, the cubic scalar vertexes among the light Higgs par- 
ticles h, H,A and h- are set equal to -z'3m^/vEw, and the cubic scalar vertexes between So- and 
two light Higgs particles are assumed to be -Z/li ^Vo-. The vertexes of //(/?, H,A) are related with 
the light Higgs mixing and the Yukawa scale factors Rq. controls the Yukawa couplings and its 
definition can be found in Eq. (28) of Ref. [19]. Once the light Higgs mixing and R^ are fixed, 
one can predict the coupling Ai d from the DM relic density. 



In the SSDM-2HBDM, the DM-nucleon elastic scattering cross section is given by [|l9h 

An m]j\mi, + mn 



-f«^/4f^^ff4+4+4f ■ (8) 
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The parameters /i , and /s have been given in Ref . [11911 and are related with the light Higgs 
mixing and the Yukawa scale factors Rq. Neglecting possible cancelation due to the light Higgs 
mixing in Eq. ([8]), we find that a^} can be enhanced by the large Rq and approach the current 
experimental upper bound for the heavy DM mass [lol. On the other hand, one can also adjust 
Rq to avoid the current direct detection limits for the light DM mass. For illustration, we take 
the Yukawa scale factors Rq = 1 for quarks and i?/ = 10 for charged leptons. Meanwhile, we 



consider the case II for the light Higgs mixing [|l9ll . In terms of the observed DM relic density, 



we calculate the allowed coupling Ai^d for 1 GeV < mo < 200 GeV. As shown in Fig. [T] (right 
panel), the observed DM abundance requires A ~ 0(10''^ - 1). Then we plot the predicted crfj^ with 



0.26 < / < 0.73 in Fig. [2] (bottom row). The latest XENONIOO [|30|] may exclude 7.5 GeV < 
ifiD < 52 GeV and 67 GeV < < 72 GeV. It is clear that the SSDM-2HBDM has smaller erf 
than that in the SSDM-SM for 1 GeV < mo < 10 GeV. In the next section, we shall see another 
advantage of the SSDM-2HBDM. Namely, the SSDM-2HBDM can give larger neutrino induced 
upgoing muon fluxes than those in the SSDM-SM even if the SSDM-2HBDM has the smaller erf. 
This is because that the large Ri can significantly change the branching ratios of the dominant DM 
annihilation channels which are relevant to the produced neutrino fluxes. 

III. NEUTRINO SIGNALS FROM THE DARK MATTER ANNIHILATION IN THE SUN 

Based on the DM mass discussed in this paper, two DM particles may annihilate into 
fermion pairs, gauge boson pairs and Higgs pairs. Therefore the diff"erential muon neutrino energy 
spectrum at the surface of the Earth from per DM pair annihilation in the Sun can be written as 

where fs denotes the DM annihilation final state and Bf^ is the branching ratio into the final state 
fs. Bfs can be exactly calculated when the couplings A and A\^d are obtained from the DM relic 
density. dN{^ /dEy^ is the energy distribution of neutrinos at the surface of the Earth produced by 
the final state fs through hadronization and decay processes in the core of the Sun. It should be 
mentioned that some produced particles, such as B mesons and muons, can lose a part of energy 
or the total energy before they decay due to their interactions in the Sun. In addition, we should 
consider the neutrino interactions in the Sun and neutrino oscillations. In this paper, we use the 



program package WimpSim [|34ll to calculate dN{^/dEy^ with the help of Pythia Il35ll . Nusigma 

10 



13611 and DarkSUSY BTll . The Pythia can help us to simulate the hadronization and decay of 
the annihilation products and collect the produced neutrinos and antineutrinos. The Nusigma is 
a neutrino-nucleon scattering Monte Carlo package for neutrino interactions on the way out of 
the Sun. The density profile of the Sun may affect neutrino oscillations due to matter effects. 
For the solar density, the WimpSim uses the standard solar model BS05(OP) [38] which is coded 
into the DarkSUSY. Notice that the WimpSim does not simulate the Higgs annihilation channel. 
Since the Higgs decay branching ratios and the energy distribution of the Higgs decay products 
can be exactly calculated in the SSDM-SM, the differential neutrino energy spectrum from the 
Higgs annihilation channel can be generated by those from other annihilation channels. Except 
for the DM masses and annihilation channels, the WimpSim only requires inputs of the neutrino 
oscillation parameters. Here we consider the lastest Daya Bay results and take [401 

sin^ 012 = 0.32, sin^ 6*23 = 0.49, sin^ ^13 = 0.026, 6 = O.S3n, 

Amli = 7.62 X 10"^eV^ Am^j = 2.53 x 10"^eV^ (10) 

for the neutrino oscillation parameters. Once dNyJdEy^ is obtained, we can use the following 
equation to calculate the differential muon neutrino flux from the solar DM annihilation: 

d^y, Tann dN,^. 



(11) 



where Res = 1.496 x 10^^ cm is the Earth-Sun distance. The solar DM annihilation rate Fann will 
be given in Eq. (1201) . In addition, we should also calculate the differential muon anti-neutrino flux 
which can be evaluated by an equation similar to Eq. ([TT). 



A. Dark matter capture rate and annihilation rate in the Sun 



The halo DM particles can be captured by the Sun via elastic scattering off solar nuclei. On the 
other hand, the DM annihilation in the Sun depletes the DM population. The evolution of the DM 
number A'^ in the Sun is given by the following equation [|4ll] : 



N = Cq-CeN- CaN^ , 



(12) 



where the dot denotes differentiation with respect to time. The solar capture rate C© may be 
approximately written as [[ill 
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where cr^ is the spin-independent cross section of the DM elastic scattering off nucleus N, . For the 
local DM density pdm and the local DM root-mean-square velocity v, we take pdm = 0.3 GeV/cm^ 
and V = 270km/s. / and 0, describe the mass fraction and the distribution of the element / in the 
Sun, respectively. /J, 0, and the form-factor suppression Ffynu) can be found in Ref. yj]. The 
function S{x) denotes the kinematic suppression and is given by 



S{x) = 

with 



1.5 l2/3 



1 +A(jc)i-5 



(14) 



3X /(Vesc)'^ 



where (Vesc) = 1156 km s~' is a mean escape velocity. In Eq. (fT2]) . the term CeN describes the 



DM evaporation rate. For the parameter Ce, we adopt the following approximate formula [|6l WM 



Ce « iQ-^-^^^-olo.vy^^-^ . (16) 

5 X lO-^^cm^ ^ ^ 

The last term CaN'^ in Eq. (fT2)) controls the DM annihilation rate in the Sun. The coefficient Ca 
depends on the thermal- average of the annihilation cross section times the relative velocity (crv) 
and the DM distribution in the Sun. To a good approximation, 

C. = <^ , (17) 



where Veit is the effective volume of the core of the Sun and is given by [14 ill 

,n ,/lGeV\^^^ 

Veff = 5.8 X 10^° cm^ . (18) 

\ rriD I 

It is worthwhile to stress that (crv) in Eq. (flTl) should be evaluated at the solar central temperature 
T, = 1.4x lO^K. 

InRefs. andH, we have calculated the DM-nucleon elastic scattering cross section cr^ 
which is equal to a^. The relation between cr^ and cr^ can be written as 

<=Al^oi;. (19) 



where An, is the mass number of the nucleus N, and M{x) = mDm^l{mo + m^). If » m^., we 
can easily derive ~ cr^. Then one may find that the solar capture rate by other elements 
in the Sun is much larger than that by the hydrogen element although it has the maximal mass 
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fraction. In terms of relevant formulas in Refs. [|l6ll and [|l9h . we calculate (crv) at = 1.4 x 10^ 
K. Using (T^ and (crv), one can straightly calculat e Cp , Ce and Ca- Then we solve the evolution 
equation and derive the solar DM annihilation rate [14 ill 



ANN 



tanh(/cro VCqCa) 



(20) 



^/^ + C£/(2 VCqCa) tanh(/cro yfC^) 

where /<• = -^1 + C|/(4CoCa) and - 4.5 Gyr is the age of the solar system. When is small 
enough (mo > 4 GeV), one may neglect the evaporation effect and obtain 



1 , I 

Tann = 2<^otanh (?o VCqCa) 



(21) 



If tQ VCqCa ^ 1, the DM annihilation rate reaches equilibrium with the DM capture rate. In this 
case, we derive the maximal DM annihilation rate Fann = Co/2 which is entirely determined by 
Cq. Therefore the enhanced (crv) via the Breit-Wigner resonance enhancement mechanism H3n 
can not affect Fann- For nio > 4 GeV, we find that most parts of the parameter space reach or 
approach the equilibrium except for the resonance region. It is because that both crf^ and (crv) are 
very small in this region [|l6ll . 



B. Neutrino induced upgoing muon fluxes in the Super-Kamiokande 

The high energy muon neutrinos from the solar DM annihilation interact with the Earth rock 
to produce the upgoing muon flux which can be detected by the SK detector [fiol . The neutrino 
induced muon flux is give by rt44ll 

ihr 1^ "/i C/j a=p,n f 

+(v^ ^ Vp), (22) 

where Pp « 1 /2NaP and p„ « 1 /2NaP are the number densities of protons and neutrons near 
the detector, respectively. Na is the Avogadro's number and p is the density of the rock under 
the detector. E^^ = 1.6 GeV is the threshold energy of the SK detector. g{L,E^,E'^)dE^ is the 
probability that a muon of initial energy E'^ has energy between E^ and E^ + dE^ after propagating 
a distance L in the rock. For the charged-current interaction cross sections, we use [4] 



d(T%{,E,,E') 2mMl( E:^^ 



dE'^ n 
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X X p 2 

) 



(23) 



whereA";" = 0.25,0.15,5"f = 0.06, 0.04 and A";f = B^y" , B"f = A'^'" . The probability ^(L,^^,^;) 
can be obtained from the full Monte Carlo calculation of muon propagation. Here we use the 



approximation formula [|44ll 



g(L,E„E') = ^^^—^, (24) 

with 

1 a+/^E' 

where a = 2.3 x 10"^ g"^ GeVcm^ and /3 = 4.4 x 10"^ g"^ cm^ describe muon energy loss in the 



standard rock H5n . It is shown that this analytic approximation is good to within 10% or better 
II44I1 . Then one can derive 

1 r™" d^y,, r^>'. ^ dcr^y(Ey^,E') 



dE, , ^ , d^^.^ dE; y ' ' Pa + (V, ^ V,) . (26) 
J^sK 'p{a+l3E^) Je^ 'dEy^JE, dE'^ 

Using a change of variable, we find that the formula in Eq. (1261 ) is consistent with that in Ref. |17|]. 

For the SSDM-SM, we calculate the neutrino induced upgoing muon fluxes in the Super- 
Kamiokande with the help of Eqs. ( fTTl) . ( [20|) and ( f26l) . The numerical results have been shown in 
Fig- in (top row). Due to the multiple Coulomb scattering of muons on route to the detector, the 
final directions of muons are spread. For 10 GeV < mo < 200 GeV, the cone half-angles range 



from 5° to 25° Il46(| . Therefore we conservatively take 0„ < 1.6 x 10 '"^cm ^s ' (maximal value 



in Fig. 8 of Ref. nl(J\ ) for the Super- Kamiokande limit. It is clear that our results in the region 
3.7 GeV < mo < 4.2 GeV and / > 0.65 slightly exceed the Super- Kamiokande limit. Since the un- 
certainties in the astrophysics and particle physics, such as puu, v and a, we can not claim that the 
Super-Kamiokande can exclude this region. Notice that the exceeded region is not consistent with 
the CDMS (shallow-site data) results as shown in Fig. |2](top left panel). For 10 GeV < mo < 200 
GeV, our numerical results in Fig. |4](top right panel) show that the predicted muon fluxes are less 
than the Super-Kamiokande limit. In Fig. 11 of Ref. ifioll . the Super-Kamiokande collaboration 
has also given the neutrino induced upgoing muon flux limits as a function of the DM mass. Their 
simulations assume that 80% of the annihilation products are from Z?^, 10% from cc and 10% from 
T^T . It is found that O^, < 6.4 x 10 '^cm^^s^' at mo = 200 GeV [lO]]. In this case, our numerical 
results are still far less than 6.4 x 10"'^cm"^s~^ 

For the SSDM-2HBDM, the large Yukawa scale factors Ri = 10 for charged leptons can signif- 
icantly enhance the branching ratio of the r^r" annihilation channel when two DM particles can 
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m„ (GeV) 



FIG. 4: The predicted neutrino induced upgoing muon fluxes for 1 GeV < mo < 200 GeV in the SSDM- 
SM (top row) and SSDM-2HBDM (bottom row). The dashed line denotes the Super-Kamiokande muon 
flux limit. 

not annihilate into PF^iy^ {nio < m^/). Since the produced muon event numbers from a pair of 
T^T" are far larger than those from bb and cc. Therefore the SSDM-2HBDM with the enhanced 
T*T~ branching ratio (5^+^- - 53% at = 10 GeV) can give larger neutrino induced upgoing 
muon fluxes than those in the SSDM-SM even if the SSDM-2HBDM has the smaller as shown 
in Figs. |2]and|4l If Ri » 10, one will obtain a smaller A\ d from the desired DM relic density 
which leads to a smaller a^. In this case, the SSDM-2HBDM will produce smaller muon fluxes 
since Ri » 10 does not significantly enlarge B^+t . For 10 GeV < mo < 200 GeV, the predicted 
muon fluxes in the SSDM-2HBDM are far less than the Super-Kamiokande limit as shown in Fig. 
m (bottom right panel). Since cr^^ in the SSDM-2HBDM may approach the current experimental 
upper bound through adjusting Rq, we can roughly evaluate the maximal muon fluxes from Figs. 
[2] and m (bottom right panels). We find that the maximal neutrino induced upgoing muon fluxes in 
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the SSDM-2HBDM are still less than the Super-Kamiokande limit when mo > mw- 



C. Neutrino induced upgoing muon event rates in tlie IceCube 

The neutrino induced upgoing muons can also be detected by the neutrino telescope IceCube 



llllh . In this subsection, we use the following formula to calculate the neutrino induced upgoing 



muon event rates in the IceCube: 

J El 2 pia+l3E^)JE^ 'dEy^jE, dE'^ 

+(v^ ^ y,), (27) 

where AesiE^) and E^^ = 50 GeV are the effective area and the threshold energy of the IceCube 
detector. To a good approximation, Agff (fi^) has a very simple functional form Hln 



A,s(E^ < lO'-^GeV) = 0, 
Aeff(10'-^GeV < < lO^-^GeV) = 0.748 [log(£'^ /GeV) - 1.6] 

Aeff(£p > lO^-^GeV) = 0.9 + 0.54[log(£^/GeV) - 2.8] kml (28) 

i?(cos 9y) is a phenomenological angular dependence of the effective area for upgoing muons 

7?(cos a) = 0.92 - 0.45 cos 9, , (29) 

where 9. is the zenith angle. Considering the change of the Sun direction, we average i?(cos 9^) 
from cos(90°) to cos(l 13.43°) and derive <i?(cos 9,)) = 1 .01 . The factor of 1 /2 in Eq. ([27]) accounts 
for about 50% of the time that the Sun is below the horizon. For the ice, we take a = 2.7 x 
10-3 g-i GeVcm2,;8 = 3.3 x lO^^'g-^ cm^p^ « 5/9A^aP and p„ « 4/9A^aP 1451]. 

We use the above formulas to calculate the muon neutrino and muon anti-neutrino induced 
upgoing muon event rates as well as the background from atmosphere neutrinos in the IceCube. 



The atmosphere neutrino fluxes d(^y^/dEy^(cos 9S) can be found in Ref. ll48h . For the atmosphere 



background, {R{co^9J)d<^v^ldEy^ inEq. (|27]) should be replaced by {R(cos9S)d<^y^/dEy^(cos9-)). 
In order to reduce the background from atmosphere neutrinos , we require E^^ < E^ < 200 GeV 
and only consider the fluxes observed along the line of sight to the Sun within the 2° half-angle 
cone isl]. Our numerical results have been shown in Fig. |5] It is found that the predicted muon 
event rates in the SSDM-SM and SSDM-2HBDM are less than the atmosphere background 10.2 
yr"^ 
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FIG. 5: The predicted muon event rates for 10 GeV < mo < 200 GeV in the SSDM-SM (left panel) and 
SSDM-2HBDM (right panel). The dashed line denotes the atmosphere background in the IceCube. 

IV. DISCUSSIONS AND CONCLUSIONS 



In terms of the observed DM abundance, we can derive the DM-Higgs couplings A in the 
SSDM-SM and ^1,0 in the SSDM-2HBDM. If and A] are enlarged by X times, the spin- 
independent DM-nucleon elastic scattering cross section crf^ in the SSDM-SM and SSDM- 
2HBDM will be enlarged by the same times as shown in Eqs. ^ and ([8]). Since the DM relic 
density will be approximately suppressed by X times, one thus needs to introduce new DM can- 
didates. In terms of Eq. (fT3l) . one may find that the produced neutrino signals from the DM 
candidates S and 5d do not significantly change as the couplings A and Ai d increase. 

In conclusion, we have investigated the singlet scalar dark matter from direct detections and 
high energy neutrino signals via the solar DM annihilation in the SSDM-SM and SSDM-2HBDM. 
Firstly, we consider the uncertainties in the hadronic matrix elements and recalculate the spin- 
independent DM-nucleon elastic scattering cross section cr,^^. It is found that the current DM 
direct detection experiments can exclude the / > 0.63 region for 1 GeV < mo < 10 GeV in the 
SSDM-SM. The latest XENONIOO may exclude 7 GeV < mo < 52 GeV and a narrow region 
65 GeV < mo < 80 GeV for m/, = 125 GeV even if we take / = 0.26. For the SSDM-2HBDM, 
we can adjust the Yukawa couplings to avoid the direct detection limits. Then we numerically 
calculate the neutrino fluxes from the DM annihilation in the Sun and the neutrino induced upgoing 
muon fluxes in the Super-Kamiokande and IceCube. The predicted muon fluxes in the region 
3.7 GeV < niD < 4.2 GeV and / > 0.65 slightly exceed the Super-Kamiokande limit in the 
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SSDM-SM. However, this exceeded region can be excluded by the CDMS (shallow-site data). We 
find that the SSDM-2HBDM can give larger muon fluxes than those in the SSDM-SM even if the 
SSDM-2HBDM has smaller a^. This is because that the large Yukawa scale factors Ri = 10 for 
charged leptons can significantly enhance the branching ratio of the t^t" annihilation channel and 
the produced muon event numbers from a pair of t^t^ are far larger than those from bb and cc. For 
the allowed parameter space of the SSDM-SM and SSDM-2HBDM, the produced muon fluxes in 
the Super-Kamiokande and muon event rates in the IceCube are less than the experiment upper 
bound and atmosphere background, respectively. The large muon fluxes in 3 GeV < nio < 10 
GeV indicate that the future neutrino experiments can provide constraints on the SSDM-SM and 
SSDM-2HBDM. 
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